Topological phases and fractional excitations of exciton condensate in a special class 

of bilayer systems 
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We study the exciton condensate in zero temperature limit in a special class of electron-hole 
bilayer systems adjacent to insulating ferromagnetic films. With the self-consistent mean-field ap- 
proximation, we find that the Rashba spin-orbit interaction in the electron and hole layers can 
induce the p ± ip or p pairing states depending on the different magnetization of the overlapped 
ferromagnetic films. Correspondingly, the topologically nontrivial or trivial phases emerge. Fur- 
thermore, in the topologically nontrivial phase, the quasiparticle excitations of the U(l) vortex are 
attached to fractional quantum numbers and obey Abelian statistics. 
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I. INTRODUCTION 



Most phases in condensed matter physics can be well 
understood with Landau's phase transition theory, which 
characterizes states of matter in terms of local order pa- 
rameters associated with spontaneous breaking of the un- 
derlying symmetries. However, the quantum Hall (QH) 
state [l], 0] gives the first example of topological order 
[H, Q| without breaking any symmetries and can be char- 
acterized by topological quantum numbers and gapless 
edge states. From then on, search for possible realiza- 
tion of nontrivial topological phases has become an in- 
triguing and challenging issue in basic condensed matter 
physics as well as technological applications. Recently, 
the two dimensional (2D) topological band insulators 
(TBI) [__-[8| present new topological phases which have 
QH effect or quantum spin Hall (QSH) effect depending 
on time-reversal symmetry (TRS) broken or not. Corre- 
spondingly, The 2D TBI can be well characterized by the 
(Thouless, Kohmoto, Nightingale, and Nijs) TKNN num- 
ber 0] for the QH phases and by Z2 topological number 
Q for the QSH phases. Soon after, the 3D TBIs, which 
are the natural generalization from 2D TBIs, were dis- 
covered in several real materials such as Bii_. E Sb :E alloys 
and Be2Se3-family crystals 

Besides 2D and 3D TBIs, more recently, a new class 
of topological superconductors (TSs) has been predicted 
by the topological classification of the Bogoliubov-de- 
Gennes (BdG) Hamiltonians [3, EH ■ Most interestingly, 
the topological quasiparticle excitations in TSs are Majo- 
rana Fermions with non- Abelian statistics [ljl [l?} , which 
have potential applications in topological quantum com- 
putation (l8j |. In analogy with TS, recently, we gave 
a proposal to realize the topological exciton condensate 
(TEC) in a spin-orbit coupled electron-hole bilayer sys- 
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tern adjacent to two insulating ferroma gne tic (FM) films 
[l9l | . Note that in our previous proposal [l9j|. the electron- 
hole bilayer is fabricated with semiconductors such as 
GaAs or InAs heterostructures, where the electron and 
hole have different effective mass and different type of 
Rashba spin-orbit coupling (namely, fc-type in electron 
layer and fc 3 -type in hole layer). Unfortunately, these dif- 
ferences break the particle-hole s ymm etry (PHS) of the 
bilayer system considered in Ref . [19| • We know that the 
PHS is indispensable for emergence of the stably topolog- 
ical quasiparticle excitations, such as Majorana fermion 
in TSs [lo, _2l| • Hence, we can predict no stable topo- 
logical quas iparticle excitations exist in TEC proposed 
111 Ref. [m. In this paper, we extend the work of Ref. 
[lH by considering the presence of PHS and discuss the 
topological quasiparticle excitations. 

We start with an electron-hole bilayer system shown 
in Fig. 1(a). In order to preserve the PHS, it is not 
suitable for the electron and hole layers to be fabricated 
with semiconductor heterostructures due to the afore- 
mentioned reasons. We argue that the recently predicted 
topological crystalline insulator (TCI) [22| films are ap- 
propriate candidates which have metallic surface states 
with quadratic dispersion on high-symmetry crystal sur- 
faces. The robustness of the surface states of TCI are 
protected by the crystal point group symmetries [for ex- 
ample four- fold (C4) or six- fold (Cg) rotational symme- 
try] associated with TRS of spinless fermion. Hence, 
the TCIs are the counterpart of Z2 topological insulators 
without spin-orbit coupling. Correspondingly, when TRS 
of spinless fermion is broken, the TCIs are the counter- 
part of TKNN topological insulators. When the suitable 
parameters are chosen in the TCI model, the dispersion of 
the surface states can be described as efe = ±k 2 /2m e ff. 
Then the electron and hole layers can be obtained by 
independently modulating the gate voltages attached to 
the TCI films. Generally, only the orbital degrees of free- 
dom are considered in TCIs. However, the spin degrees 
of freedom have to be considered if the magnetic effect 
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FIG. 1: (Color online) (a) Schematic structure of the effec- 
tive electron-hole bilayer system that holds the exciton con- 
densate. Different layers are fabricated with corresponding 
material shown on the layers, in which the spacer are fab- 
ricated with dielectric and its thickness is d. The electron 
and hole layers are obtained by tuning the gate voltages V ge 
and V g h- The inserted picture is the dispersion of the sur- 
face state of TCI. The FM films provide effective exchange 
fields (V e ,Vh)- (b) The energy spectrum of the effective 
electron-hole bilayer system near the Fermi energy Ef=0. 
Here the dashed lines denote non-interacting single-particle 



energy spectrum E t 



(h)±, 



while the solid lines denote the exci- 



ton energy spectrum with an obvious mean-field gap opened. 
We take t e = t h = 1, Ate = iih = -4, (V e , V h ) = (-1, -1), and 
a = 0.5. 



is introduced, which is the case in our model. We intro- 
duce two insulating FM films adjacent to the two TCI 
films to support two effective interface exchange fields, 
which respectively induce the population imbalance of 
the spin-up and spin-down carriers in the electron and 
hole layers and break the spinful TRS of the system due 
to the Zeeman effect of the exchange fields. We argue 
that the magnetic vector potential (orbital effect) of the 
exchange fields of the FM films can be neglected, be- 
cause the exchange fields only locate on the interface of 
the FM films and TCI films. That is the reason why 
we do not apply external magnetic fields to replace the 
FM films. Additionally, we assume that the exchange 
fields are always perpendicular to the bilayer. It is pre- 
dictable that the interlayer Coulomb interaction induces 
the excitonic instability when the two layers are sepa- 
rated with an insulating spacer and close enough with 
each other. Moreover, We note that since structure in- 
version symmetry is broken at the TCI film and the insu- 
lating spacer interface, there will be an interface Rashba 
spin-orbit coupling (SOC) [H SI]: U R = a (a x k) • z 
and a cx (V Z V) , where V comes from doping, impuri- 
ties, external gate voltages and so on. It is evident to find 
that a is not zero when the structure inversion symmetry 
is broken and can be tunable by external gate voltages. 
Moreover, a Dresselhaus term [25| can not be zero when 
the bulk inversion symmetry is broken. In the two dimen- 
sions, Hd — P{k x a x — k y (jy) and (3 ~ —B(ir/w) 2 , where B 
is a material-dependent parameter. In the conventional 
semiconductor quantum-well systems, w is the width of 
the well with a typical value of ~50 A, then /3~0.3 (0.1) 
eV A in InSb (GaAs) quantum wells, which is compa- 



rable to or even larger than the value of a [26(. In this 
case, the Dresselhaus term can not be neglected. How- 
ever, if the electron and hole layers are fabricated with 
TCI films, since the TCI has the bulk inversion symme- 
try, then the Dresselhaus term can be neglected and we 
can only consider the Rashba term. When the amplitude 
of the Rashba SOC is comparable to that of the exciton 
order parameter (EOP), we find that the unconventional 
pairing states emerge with topologically nontrivial p+ip 
pairing states or trivial p pairing states depending on 
the magnetization of the overlapped FM films. Further- 
more, we find that the topological quasiparticle excita- 
tions of the U(l) vortex profiles survive and carry frac- 
tional quantum numbers [13, S3 that obey the Abilan 
statistics. 

The paper is organized as follows. In Sec. II a model 
Hamiltonian for the effective electron-hole bilayer system 
is introduced, and a self-consistent mean-field method is 
used to produce the unconventional pairing states. In 
Sec. Ill we introduce the TKNN number as the topolog- 
ical number to characterize topologically nontrivial pair- 
ing states and show the chiral gapless edge states. In Sec. 
IV we analytically study the topological quasiparticle ex- 
citation of a single U(l) vortex in the low-energy limit 
and give the numerical results in lattice model and find 
both are consistent. In Sec. V we summarize our results. 



II. MODEL HAMILTONIAN AND EXCITON 
CONDENSATE 

Consider the effective TCI electron-hole bilayer con- 
fined in the x-y plane in zero temperature limit. Then 
the mean-field description is appropriate even in two di- 
mensions. For convenience of discussion, here we con- 
sider a square lattice model. Furthermore, the contin- 
uous model can be obtained from the low-energy limit 
of the lattice model. The square lattice Hamiltonian is 
u—^r cub) ^u^ pS >\^T-i^ e ~ h ^ ■ 



H 



(p) 

kin 



tlr, — 



nj(e-h) 
rL int 



(^P ~ to^PtcrPj,* 

<i,j>,u 
j 

3 

- U i , j (d)el <T h t ja ,h jtr >e itT . (1) 



Here t p denotes the nearest-neighbor hopping amplitude 
while /i p represents the chemical potential in electron 
(p=e) or hole (p=h) layer. Without loss of generality, 
we set t p =l for both electron and hole layer, other pa- 



rameters are measured in t„ 



are the three Pauli 



3 



matrices and sq is the identity matrix. V p represents Zee- 
man splitting from the effective interface exchange field 
of the FM film. pj t(T is the annihilation operator of car- 



rier with spin a at lattice site j. 



is the Rashba 



SO interaction strength in the electron (hole) layer. For 
simplicity, in this paper we set a e =ah=a. Sx (Sy) is the 
square-lattice spacing measured in unit of in-plane lat- 
tice constant along the x (y) direction. The interlayer 



Coulomb interaction is Uij(d)=e 2 /ey\ri te — r^j + d 2 , 
where e is the dielectric constant of the spacer and d is 
the interlayer distance. We only consider the interaction 
relevant to exciton formation and ignore the electron- 
hole exchange interaction. We also neglect the intralayer 
electron-electron and hole- hole interactions, since they 
are expected to renormalize the single-particle spectrum 
of each layer and have no essential influence on the topo- 
logical properties of the system. 

In momentum space, the Hamiltonian Eq. (1) can be 
expressed as 

= y^a(sin k x - i sin fc^pj^Pk; + H.c, 



ri int 



E 

kk'qtrtr' 



q k+qrr k' — q<r 



^k'-q<T'^k' CT 'e kCT , (2) 



where U(q)= 2 ^-e-'' d and Q p, =-2i p (cosfc a; +cosfc a ). In 
the Nambu notation with combined e-h field operator 
basis ^k=[ekt ^kj. ^l k t ^-kJ T > ^ ne decoupled mean-field 
Hamiltonian is expressed as Hmf=^'H'Ho + 'Hi)^+Eq in 
the following Bogoliubov-de Gennes (BdG) form: 



-(e) 



- fX e + V e S z 

At(k) 



J -k 



A(k) 

+ Mfc - VhS z 



( 3 ) 

where T~Lq represents the single-particle terms while Tii 
includes the interaction terms, E^=±£2k s o±%if' > (±k), 

A ^'( k ) = -sE q t / (q)(' l -k+q ff 'ek- qff ) are the EOPs 

and E Q - 



lv A_,(k)A^,(k- q ; 



in which we have 



neglected the inessential constants. 

Moreover, the exciton condensate with unconventional 
pairing can be well understood in the Fermi surface 
nesting picture. To reveal this fact, firstly, the non- 
interacting Hamiltonian "Ho i n Eq. Q can be diagonal- 
ized with the unitary transformation in the single-particle 
eigenstate space as 



= E*Sk (#i e) (k) © 4 h) (k)) *ok, (4) 

k 

where H$ (k)=E p _ (k)(s + s z )/2+E p+ (k)(s Q -s z )/2, 
* ok = [4+ (k) , 4- (k) , ft h+ (k) , $ _ (k)] T , 



S ps (A;)=±c£ 3 2=FMp+s\/ a 2 (sin 2 sin 2 
(s=+, — ) are respectively electron and hole band ener- 
gies, and ipps denote the relevant particle annihilation 
operators. The two basis sets are related by the corre- 
sponding unitary transformation, V& \=e l9k Ii e © e k n^ 
^ok with 



n„ 



where 



-if p+ (k)e-^ f p .(k) 
/ p _(fc) -if p+ (k)e^ 

^> k =arctan(sin k y / sin k x ) 



/p±(k)=Wk/ 



v 2 ± v p y 



(5) 

and 
with 



Lo-^—a^ sm 2 k x + sin 2 k y . It should be stressed that 

the fc-dependent phases 9k and are confirmed by 
exactly solving the ground state of the system through 
self-consistent mean-field calculation. The single-particle 
bands E p ±(k) are shown in Fig. 1(b) (dashed curves), 
in which the occupied bands are E e _ and Eh+. With 
the values of the tunable parameters shown in the 
caption in Fig. 1, the band E e _ and band Eh+ have 
the perfect identical Fermi surface. From the nesting 
mechanism, the attractive interlayer interaction leads 
to unstability of the Fermi surfaces and opens a gap. 
This is the case similar to the conventional BCS picture. 
Now, the electron-hole interaction part relevant to 
exciton formation in Eq. ^ can be expressed in terms 
of the occupied electron band E e ^ and hole band 
Eh+ with Eq. §5§. Then, the decoupled mean-field 
two-band Hamiltonian of the exciton system under the 



basis \&= 



constants neglected) 



V> e _(k), tM. + (— k) reads (with inessential 



MF 



E* k 



E e - A(k) 
A*(k) E h+ 



Eq, 



(6) 



where A(k)=-£ E U ™0*> s ' s ') (V^+(-k+q)?A e - (k-q)) 

with renormalized in- 



q.s,.s*' 
A(k)A*(k') 



and £ =i^ 

kq 2^ !7 *«( k ' C l' S ' S ' > 

teraction U re (k, q, s, s')=U (q)xk, q ^ SS ' (k, q)r k , q , 

s=±, Xk]q = e i(ek-<,-flk) e i(i?k-^k- q ) 5 an d 

F ss >(k,q)=f es (k)f es (k-q)f es ,(k)f es ,(k-q). Here, 
Tk, q =e i(0k - C1 ~ 0k) for s^s', r k , q =l for s=s'=-, and 
7- k , q =e i2(0k - C|_0k) for s=s'=+. Then the EOPs in Eq. 
(O and those in Eq. ^ are related by the unitary 
transformation associated with Eq. ([5]), i.e., 



A(k) = -e l(ek - ,?k) A(k) 

if e+ {k)f h+ e-^ / e+ (k)/ A _(k) 
-/ e _(k)/ h+ (k) i/ e _(k)/ h _(k)e^ 



(7) 



The factors / e _|- (k)/^± (k) in Eq. ([7} can be understood 
as follows. Since the four pairing states with different 
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FIG. 2: (Color online) The self-consistent results of EOPs ver- 
sus momentum (k x , k y ). (a), (b), (c) and (d) refer to |A^(k)|, 
|A n (k)|, |A4. t (k)| and A u (k)| .Correspondingly, (al), (bl), 
(cl) and (dl) refer to Re(A n (k)), Re(A tt (k)), Re(A it (k)) 
andRe(A u (k)); (a2), (b2), (c2) and (d2) refer to Im(A n (k)), 
Im(A n (k)), Im(A it (k)) and Im(A u (k)). The effective Zee- 
man fields are (V e , Vh) = { — 1, — !)■ 
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FIG. 3: (Color online) The coefficients / e+ (k)/ h+ (k), 
/ e+ (k)/ h _(k), / e _(k)/ h+ (k) and / e _(k)/fc_(k) versus 
the momentum k in the magnetization configuration 
(V e ,Vh)=(-l, -1} for different Rashba SOC values: (a) 
a=0.5; (b) a=2. 



spin combinations in EOPs are affected differently by the 
effective exchange fields. Additionally, the Rashba SOC 
flaws the spin polarization of the carries along the z direc- 
tion, then the total effect leads the factors / e ±(k)//i±(k) 
to emerge in four pairing states of EOPs and deter- 
mine the EOP behavior at given (V e , Vh,)- We find that 
only one component of A CT(T / (k) is dominant, while the 
other components can be neglected in the proper range 
of the Rashba SOC parameter for the given magnetiza- 
tion of the FM films. Moreover, the system's proper- 
ties are decided by the dominant component of A CTCT ' (k). 
For definiteness, we have carried out a series of self- 
consistent numerical solutions of the BdG equation (3) 
by setting the lattice size 81x81, fi e =fih=—4, a=0.5, and 



«-*x)=1.0. In the 
which means the 



(V e , V/t)=(±l, ±1). As a typical example, we show in Fig. 
2 the numerical results of A aa > (k) for (V e , Vft)=(— 1, —1). 
Also, we present in Fig. 3(a) the results of the factors 
/e±(k)/ft±(k) for comparison. In this case, from Fig. 2 
it is clear that A^(k) is the dominant EOP component 
which has the phase of — ie _l ^ k , i.e. e 1 ^ 
continuous limit, — ie~ l * k ~ — "t. 1 * , 
p-Hp-like pairing state. Furthermore, this kind of pairing 
state is usually topologically nontrivial (see discussion be- 
low). Similarly, for the case of (V e , V/ l )=(l, 1), the Aj,j,(k) 
turns out to be the dominant pairing component with a 
phase of —ie l ^ k , which thus is also topologically nontriv- 
ial. On the other hand, for the case of (V e , Vh)=(— 1, 1), 
the A^(k) is the dominant pairing component with a 
phase of — (e l * k +e~ I ^ k ). The physics related to this kind 
of phase is that e l ( e k--&w) — e ±i<Pk correspond to two de- 
generate states and their linear combinations give the 
ground state with the phase of — (e l ^ k +e~ l ^ k ). In the 



continuum limit, — (e 



) 



k ■ 



This kind of 



pairing state is topologically trivial. Similarly, for the 
case of (V e ,Vh)=(l, —1), the dominant pairing compo- 
nent A^(k) is also topologically trivial. 

Furthermore, from Fig. 3 (a) and Eq. (7) we can 
find that the weak Rashba SOC relative to the Zeeman 
fields will select one dominating component of EOPs 
(namely, A^(k)), while other components can be ne- 
glected. But in Fig. 3 (b) we can find that the ampli- 
tudes of f e+ (k)f h+ (k), f e +{k)fh~{k) and f e -(k)f h+ (k) 
are comparable. That means the non-weak Rashba SOC 
relative to the exchange field cannot select one domi- 
nating component of EOPs, and all components of EOPs 
play important roles. The consequence is that the uncon- 
ventional pairing states become unstable in case of Fig. 
3 (b), where the value of a is chosen to be much larger 
than the exchange filed. We can also clarify it from nu- 
merical calculation. Hence, the Rashba SOC cannot be 
larger than the Zeeman splitting in our system. Fortu- 
nately, it is easy to satisfy this requirement, because both 
Rashba SOC and Zeeman fields are tunable in our sys- 
tem. Besides, another essential point we stress is that 
the topologically trivial or nontrivial EC in our system 
is stable, which means a stable bulk gap separating the 
ground state and the excited states. We can obviously 
find this fact in Fig. 1(b) and Fig. 4(a) (see below). In 
the following section, we focus on the characterization of 
the topological properties of the system. 



III. 



CHIRAL TOPOLOGICAL ORDER 



From the self-consistent results of the aforementioned 
section, we can find that the EOPs are closely cor- 
related to the magnetization configuration of the ex- 
change fields. We also note that there are two different 
classes: (V e , V h )=(±l, ±1) and {V e , V h )=(±l, =Fl)- So we 
choose one configuration from each class as an example. 
The EOPs can be simplified as A-j-f (k)= — iAo(sin k x — 
ismky) and other A CT(T >(k)^0 for (V e , V/,)=(— 1, — 1), 
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while A-|-j_(k)=— Ao sin k x and other A CT(T /(k)~0 for 
(V e , Vh)=(— 1, 1). Then the mean- field Hamiltonian Eq. 
(|3]) in these two magnetization configurations can be adi- 
abatically mapped into the following forms, which are 
topologically equivalent to initial mean-field Hamiltonian 
Eq. ([3]) when a — >• without the bulk gap closing [29J. 



MF 



-kt 



ff'(k) 



-kt 



(8) 




(b) 



'kt 



h- 



(9) 



where iJ*(k) = 
Vh)(«o - «s)/2 



(C k e) - 



■/x e + \4)(so + S3)/2+(- 



.(ft) 
_( .k " 



Mft 



iAo(sinfc :c si + sinfcj,S2) and _ff s (k) 



(c k e) 



Me + K)(so + s 3 )/2 + (~ci h) + »h + V h )(s - 
S3)/2 — Ao sinfczSi. It is easy to check that both H s (k) 
and H l (k) preserve the PHS and break the TRS accord- 
ing to the definition of Ref. [30| ■ The missing of TRS can 
be traced back to the effective Zeeman fields in Eq. ([3]). 
From Ref. [30( , we can find that the Hamiltonian describ- 
ing our system generally belongs to class D in BdG class. 
The topologically nontrivial pairing has the form of chiral 
p wave (namely, p x +ip y or p x —ip y pairing) in two dimen- 
sions. The p x or p y pairing is trivial. Furthermore, the 
nontrivial topology is characterized by an integer which 
is similar to the TKNN number in the QHE. 

We can calculate this topological number from Eq. ([3]) . 
The topological nature of the ground state of the exci- 
ton condensate can be characterized by non-zero Itknn, 
which reads 



1 



ITKNN 



2 

71 = 1 



d 2 kQ n (k), 



(10) 



BZ 



where fi n (k)=— 21m 



pM„(k) 
<9fc x 



du n (k) \ 
dk v I 



is the Berry cur- 



vature of the occupied energy band n and|u n (k)) is the 
corresponding Bloch wave function. The straightforward 
calculation gives Itknn=^ f° r 04, V|,)=(— 1, — 1) and 
Itknn=0 for (V e , Vh)=(— 1, 1). From the bulk-edge cor- 
respondence, the nontrivial bulk topological number im- 
plies gapless edge states emerging in the terminal of the 
system. 

In order to evaluate gapless edge states, we assume 
that the square-lattice system has two edges in y direc- 
tion and is boundless in x direction. Correspondingly, 
we choose open boundary condition in y direction and 
periodic boundary condition in x direction in the lat- 
tice Hamiltonian in Eq. (JXJ) in mean-field approxima- 
tion. The calculated energy spectrums of the topologi- 
cally nontrivial cases of (V e ,Vh)=(— 1,— 1) and topolog- 
ically trivial case of (V e , Vh)=(— 1, 1) are illustrated in 
Fig. 4 (a) and (b), respectively, in which the red-solid 
and blue-dashed lines correspond to the different edge 
states with opposite chiralities. It is obvious to find that 



FIG. 4: (Color online) The energy spectrum of the bilayer 
square-lattice system with two edges at the y direction. k x 
denotes the momentum in the x direction. The magneti- 
zation parameters are set at (V e , Vh)=(— 1, — 1) for (a) and 
{V e ,V h )=(-l,l) for (b). In both cases A =0.5 and a=0.5. 
The red-solid and blue-dashed lines denote the edge states 
locating at different edges. 



the number of the gapless edge states is consistent with 
the bulk theory characterized by Itknn- 

Furthermore, for the case of (V e , V/,)=(— 1, — 1), the 
edge states can be described by the effective Hamilto- 



H edge = ± ^2 v Fk x j\k x ,y)^(k x ,y), (11) 

k*>0 

where ± represents the opposite chiralities, vf is the 
Fermi velocity, and k x is the momentum measured from 
the Fermi surface, and 

l{k x ,y) = u(k x ,y)e f (y) + v(k x ,y)h\(y). (12) 

We will show in next section that the low-energy topolog- 
ical quasiparticle excitations are dominated by the edge 
states. 



IV. FRACTIONAL QUSIPARTICLE 
EXCITATIONS IN EXCITON CONDENSATE 

We now demonstrate that there exist fractional excita- 
tions in our system in the topologically nontrivial phases 
if the EOPs have a U(l) vortex profile with odd vorticity. 
Note that the even vorticity is trivial [3l|. For definite- 
ness, we consider the case: (V e ,Vh)=(—V, — V)=(— 1, — 1) 
and /i e =/z/j=/x=— At p . The non-vortex BdG Hamiltonian 
in Eq. (|3|) reads 



7~Lmf (k) = 



-(e) 



- n - Vs z 

At(k) 



A(k) 



Vs 2 



(13) 



with 
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A ,,, \ -iA (sink x - is'mky) 0] , 1yn 
A ( k ) = o ' ( ' 

The PHS of U MF {k) is defined by 
AWjvfir(k)A=— Hmf(— k) with A=s x ® s £, where 
/C is the complex-conjugation operator and A 2 =f . 

Now, we consider a vortex in the EOPs, i.e., the uni- 
form Ao in Eq. (14) is modulated to f 2 (r)e mS , where r 
and 6 are polar coordinates centered on the vortex, n is 
the vorticity, and f(r) is a real function of r that van- 
ishes at small r and obtains asymptotic value /o (a real 
constant) at large r. Usually, the nature of the excitation 
bound to the vortex is dominated by the low-energy limit 
of the system, i.e. k— H). In this limit, the effective Hamil- 
tonian can be obtained by the replacement k — > — i V and 
is accurate up to C(k 2 ). The zero- mode equation reads 
Hmf(-«V)* o =0 with * o (r,60=[w t ,u 4 .,v T ,t;. l ] T With 
the help of AHmf{— «V)A=— Hmf(— iV), it is easy to 
find v a =u%. Hence the zero-mode equations can be ex- 
pressed as 




FIG. 5: (Color online) The density distributions of quasipar- 
ticle excitation of the U(l) vortex locating on the center of a 
31 x 31 lattice for (a) u-|-(r)| 2 , (b) |tq.(r)| 2 , (c) i^(?-)| 2 , and 
(d) |uj.(r)| 2 in the exchange fields (14, Vh) = (— 1, —1). We take 
/ 2 =0.5, a=0.5. 



y dzdz* 



dz* 



-a— 1 + (-__ + %=0, 
oz ozaz* 



(15) 

<ui = 
(16) 



where z=x+iy. From Eq. (|12p in Sec. IV, we can find 
that can be neglected. Moreover, the numerical results 
shown in Fig. 5 also give a proof. Hence, we only need 
to consider Eq. (fl~5"|) . For dcfinitencss we consider the 
simplest case in the quantum limit analogous to the Lon- 
don approximation in superconductor, i.e., /(r=0)=0, 
/(r>0)=/o (a positive real constant) and the U(l) vortex 
with vorticity n=l. Moreover, we look for a solution in 
terms of a spherically symmetric real function u^. Note 
that in Eq. (TTB1) . we must take into account the U(l) 
symmetry in order to get the physical solutions. Then 
Eq. (TTSl) can be simplified into the following differential 
equation: 



(17) 



The solution of Eq. 



Jo(\/^V 2 + /of/2)e" / " r/2 where J a (x) is the Bessel 
function. We also give the exact numerical results of 
|u CT | 2 and \v a \ 2 in the square-lattice system, which are 
shown in Fig. 5. Clearly, it can be seen that u^O, 
which is consistent with Eq. (TT^|) . 

The quasiparticle excitation field associated with 
the zero-mode solution in the exciton condensate can 
be expressed as % C =J dr(itf&|-(r) + v^-ht(r)). Con- 
trastively, the quasiparticle excitation field describing 



the vortex zero mode in topological superconductor is 
jsc= J dr(ufdf(r) +VfcL(r)) [2lll33|. It is obvious to find 

due 



that % c ^ y| c due to [e t (r), h\(r)] + = 

to [c t (r),c|(r)]+=l. That means the quasiparticle exci- 
tations in topological superconductor are the Majorana 
fermions which obey the non-Abelian anyonic statistics 

17]. In our case, the quasiparticle excitations are not the 
Majorana fermions and are charged. Generally, a local- 
ized zero mode in a system with PHS is known to carry 
a fractional charge of ±e /2 [27l . l33j . Based on the argu- 
ments given in Refs. [34l |35|. we expect for the present 
system that the quasiparticle excitation bounding to vor- 
tice has fractional charge e/2. We can also evaluate it 
from the spatial integration Jrfr(|u-|-| + \vf \ ) around 
the vortex defect. Fractional charged quasiparticles con- 
fined to 2D often obey fractional exchange statistics [36[ , 
which is in our case Abelian from the standard arguments 

371. 



=0 and j sc = 7+ c 



V. CONCLUSION 

In conclusion, we have studied the exciton condensate 
in zero temperature limit in a special class of effective 
TCI bilayer system which is different from the semicon- 
ductor bilayer system in Ref. [19j due to its preservation 
of PHS. With the self-consistent mean-field method, we 
have found that the Rashba SOC in the electron and 
hole layers can induce the unconventional p ± ip or p 
pairing states that depend on the magnetization of the 
FM insulating films. Correspondingly, the p ± ip pair- 
ing states represent topologically nontrivial exciton con- 
densate phase, which can be characterized by nonzero 
TKNN number or the gapless edge states, while the p par- 



7 



ing states are topologically trivial. Furthermore, by low- 
energy analytic solution as well as exact numerical calcu- 
lation, we have shown that in the topologically nontrivial 
exciton condensate phase, the presence of U(l) vortexes 
result in the simple Fermi zero-mode fractional quasi- 
particle excitations that obey Abelian statistics, which 
completely distinguish from the non-Abelian Majorana 
Fermi zero-mode quasiparticle excitations in topological 
superconductor systems. 
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